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NOMENCLATURE 


length  of  converging-diverging  portion  of  test  section,  M.9  m 
pressure  coefficient,  2(P-Pi)/pU2i 

A 

=  x  /pU  2,  local  skin-friction  coefficient 
w  e 

frequency,  Hz 

=  f  +  f  +  f  ,  instantaneous  signal 

s  f  +  f  ,  ensemble-average  of  F  for  a  specific  phase 

A,  A 

shape  factor  =  6*/0 
thermal  conductivity 
=  toC/2Uie>  reduced  frequency 

distance  from  the  wall  to  the  minimum  velocity  in  the  ensemble 

averaged  velocity  profile  backflow 

pressure 

=  U^e/v ,  momentum  thickness  Reynolds  number 

period  of  imposed  oscillation 

time 

time-averaged  velocity 
Periodic  velocity 

ensemble-averaged  velocities  in  x,  y  directions,  respectively 
turbulent  fluctuations  in  x,  y  directions,  respectively 
ensemble-averaged  Reynolds  shear  stress 
ensemble-averaged  mean  square  values  of  u  and  v 

-  (x  /p)  ^  2 ,  shear  velocity 
s  U/U 

T 

maximum  backflow  velocity 

mean  test  section  entrance  velocity 


U  amplitude  of  n  harmonic 

ne 

x,y,z  streamwise,  normal,  and  spanwise  coordinates 
+ 

y  a  Reynolds  number  =  ylWv 

Greek  Symbols 


6 

6* 

< 

0 


rne 


pu 


pu 


min 


ie 


OJ 

Subscripts 

e 

n 

w 


y  where  J  =  0.99  U 

e 

displacement  thickness 
von  Karman  constant,  0.^1 
momentum  thickness 

phase  angle  of  nth  harmonic  of  velocity 
kinematic  viscosity 

ensemble-averaged  forward  flow  fraction  at  1  mm  from  the  wall 

Y  at  a  distance  N  from  the  wall 
pu 

phase  angle  of  first  harmonic  pressure  gradient 
density 

shearing  stress 
2irf 

freestream  or  condition  outside  shear  flow 
harmonic  number,  e.g.,  1,  2,  etc. 
denotes  wall  value 


min, max  minimum  and  maximum  value 


1 .  INTRODUCTION 


Unsteady  turbulent  boundary  layers  are  of  considerable  interest  because 
of  unsteady  aerodynamic  phenomena  associated  with  blades  in  compressors  and 
turbines,  helicopter  rotors  in  translating  motion,  high  pressure  recovery 
diffusers  and  aeroelastic  flutter  and  support  of  aerofoils.  They  are 
particularly  important  during  high  loading  conditions  when  separation  may  be 
present  during  a  portion  of  the  oscillation  cycle.  Under  such  conditions 
significant  interaction  between  the  thick  turbulent  shear  layer  and  the 
inviscid  flow  occurs.  Thus  experimental  information  on  the  structure  of 
separating  unsteady  turbulent  shear  layers  is  necessary  for  the  proper 
understanding  and  calculation  of  the  physical  behavior  of  these  practical 
flows. 

Simpson,  Shivaprasad  and  Chew  (1983a)  and  Simpson  and  Shivaprasad 

(1983b)  have  reported  the  measurements  showing  the  effect  of  frequency  on 

moderate  amplitude  (amplitude  to  mean  velocity  ratio  of  0.33)  sinusoidal 

unsteadiness  of  the  free-stream  velocity  on  a  separating  boundary  layer  at 

reduced  frequencies  k  =  wC/2U  .  of  0.61  and  0.90,  where  U  .  is  the  test 

e  i  ei 

section  entrance  time-averaged  free-stream  velocity,  C  =  4.9  m,  length  of 
the  converging-diverging  test  section  shown  in  Figure  1,  and  u  =  2nf,  where 
f  is  0.596  Hz,  the  frequency  of  flow  oscillation. 

In  this  report  measurements  are  presented  for  two  large  amplitude 
(amplitude  to  mean  velocity  ratio  of  0.8)  gust  flows  produced  by  an  in- 
stream  rotating  blade  damper  at  k  =  0.61  and  1.33  and  an  oscillating  flow 
with  k  =  1.03  produced  with  a  test  section  "roof  damper".  Detailed 
ensemble-averaged  velocity  and  turbulence  measurements  were  made  using  hot¬ 
wire  and  laser  anemometers.  While  the  primary  function  of  this  report  is  to 
present  and  document  these  data,  some  discussion  of  the  results  and  a  number 


of  data  plots  are  presented.  The  reduced  data  are  available  on  computer 
diskettes. 

Agarwal  and  Simpson  (1986,  1987a, b)  have  reported  some  results  for 
the  two  large  amplitude  flows.  These  large  amplitude  flows  show  some  unique 
features  not  observed  in  moderate  amplitude  flows,  viz.,  during  strongly 
decelerating  and  minimum  free-stream  velocity  phases,  backflow  occurs  in  the 
diffuser  section  and  even  near  the  wall  in  the  converging  section.  This 
appears  to  be  due  to  the  strong  oscillating  pressure  gradient  that  opposes 
the  slow-near-wal 1  flow  to  produce  backflow.  In  the  absence  of  flow 
reversal,  as  in  the  moderate  amplitude  flows  studied  earlier,  the  flow  is 
quasi-steady.  It  is  clear  from  earlier  results  of  Simpson  et  al .  (1983)  and 
those  for  these  flows  that  substantial  hysteresis  of  the  turbulence 
structure  occurs  for  unsteady  separating  flows  and  the  oscillation  waveform 
and  amplitude  strongly  affect  the  detached  flow  behavior.  New  information 
on  Reynolds'  shear  stress  and  phase  lag  of  the  turbulence  structure  with 
respect  to  the  ensemble-averaged  velocity  profile  structure  has  been 
obtained  and  is  discussed  briefly  here.  More  detailed  discussions  will  be 
contained  in  journal  publications  (Agarwal  and  Simpson,  1988a, b). 

2.  EXPERIMENTAL  EQUIPMENT 

2.1  Basic  Wind-Tunnel 

Figure  1  is  a  side  view  schematic  of  the  8  m  long,  0.91  m  wide  test 
section  of  the  wind  tunnel.  The  test  wall  is  constructed  from  18  mm  thick 
fin-form  plywood,  reinforced  every  28  cm  with  7.6  x  3.8  x  0.6  cm  cross- 
section  steel  channel.  This  reinforcement  was  necessary  since  Acharya  and 
Reynolds  (1975)  found  that  test  wall  amplitudes  as  small  as  0.025  mm 
produced  up  to  10?  error  in  U  measurements  near  the  wall.  However,  some 
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test-wall  oscillation  was  still  present,  which  was  removed  by  supporting  the 
test  wall  directly  under  the  point  of  measurements,  or.  columns  attached  to 
concrete  floor.  The  side  walls  of  the  wind  tunnel  are  made  of  float  glass 
to  prevent  laser  signal  dispersion,  while  the  upper  wall  is  made  o^ 
plexiglass. 

The  mainstream  flow  of  the  blown,  open-circuit  wind  tunnel  is 
introduced  into  the  test  section  after  passing  through  a  filter,  blower,  a 
fixed-setting  damper,  the  in-stream  rotating-blade  damper  discussed  in 
section  2.2  below,  a  section  of  honeycomb  to  remove  the  mean  swirl  of  the 
flow,  seven  screens  to  remove  much  of  the  turbulence  intensity,  and  finally 
through  a  two-dimensional  4:1  contraction-ratio  nozzle  to  further  reduce  the 
longitudinal  turbulence  intensity  while  accelerating  the  flow  to  the  test 
speed  (see  Figure  2). 

The  inviscid  core  flow  was  uniform  within  0.05X  in  the  spanwise 
direction  and  within  1 1t  in  the  vertical  direction.  The  test  wall  boundary 
layer  was  tripped  by  the  blunt  leading  edge  of  the  plywood  floor,  the  height 
of  the  step  from  the  wind-tunnel  contraction  to  the  test  wall  being  0.63  cm. 

To  obtain  additional  information  on  how  oscillating  pressure  gradient 
affects  the  detached  flow  structure  an  unsteady  flow  generated  with  a  damper 
in  the  roof  of  the  test  section  and  with  k  =  1.03  was  also  studied. 

In  this  flow  case  the  first  8  foot  section  remained  unchanged  from  the 
arrangement  used  for  all  other  flows.  The  roof  of  the  second  or  diverging 
section  was  made  from  a  fine  screen  for  producing  an  easily-controlled 
pressure  gradient.  A  set  of  12  rotating  blade  dampers  mounted  above  this 
section  (called  "roof  damper")  controlled  the  phase  and  amplitude  of  the 
part  of  the  main  flow  diverted  through  the  roof.  A  perforated  plate  at  the 
test  section  exit,  as  used  by  Saripalli  and  Simpson  (1980),  was  used  to 


increase  the  test  section  back  pressure.  A  wooden  piece  0.38  m  long  was 
placed  over  the  roof  to  delay  detachment.  This  arrangement  is  schematically 
shown  in  Figure  3.  This  flow  will  be  referred  to  as  k  =  1.03  flow  with 
"roof  damper".  For  this  flow  the  in-stream  rotating  damper,  used  for  the 
large  amplitude  flows  was  kept  fully  open  and  stationary. 

2.2  Programmable-Rotating-Blade  Damper 

Investigators  of  experimentally-produced  unsteady  flows  normally  have 
little  control  over  the  waveform  of  the  flow  unsteadiness.  Although  the 
waveform  may  consist  principally  of  a  given  frequency,  substantial 
contributions  normally  come  from  higher  harmonics.  The  in-stream  rotating- 
blade-damper  system  used  by  Simpson,  Shivaprasad  and  Chew  (1983)  and  Simpson 
and  Shivaprasad  (1983)  and  described  in  detail  by  Simpson,  Sallas  and 
Nasburg  (1978)  was  also  used  here  for  the  two  large  amplitude  (k  =  0.61  and 
1.33)  gust  flows.  In  this  system,  the  angular  velocity  of  the  rotating 
blades  in  the  in-stream  damper  in  the  flow  downstream  of  the  blower  is 
varied  during  a  cycle  so  as  to  produce  the  desired  waveform  shape,  amplitude 
and  frequency  of  the  flow  entering  the  test  section.  In  the  current 
experiments,  the  0.596  Hz  and  0.95^  Hz  flow  oscillations  had  a  velocity 
amplitude  of  about  ~  0.8  of  the  mean  velocity. 

The  0.95^  Hz  flow  oscillation  produced  with  the  "roof  damper"  had  a 
velocity  amplitude  of  ~  0.12  of  the  mean  velocity.  In  this  case,  the  in- 
stream  damper  in  the  flow  downstream  of  the  blower  was  not  used. 

The  electronics  were  the  same  in  all  experiments.  All  events  during  an 
oscillation  cycle  were  synchronized  with  respect  to  a  "reference"  square 
wave  voltage  signal  at  the  oscillation  frequency  generated  by  the  quartz 
clock  in  the  control  electronics.  A  'clock'  square-wave  voltage  signal  with 


a  frequency  96  times  the  reference  signal  is  also  generated  to  aid  data 
acquisition.  Data  were  acquired  at  the  beginning  of  each  of  96  bins  of  a 
cycle.  Reference  signal  was  also  recorded  to  ensure  the  phase  repeatability 
of  the  acquired  data. 

The  variation  in  the  period  of  each  flow  cycle  and  the  'jitter'  follows 
a  Gaussian  distribution  and  Simpson  et  al .  (1983)  estimated  that  jitter  has 
no  affect  on  ensemble-averaged  velocity  values  and  only  an  effect  on  low 
turbulence  level  fluctuation  measurements. 

2.3  Boundary  Layer  Control  System 

The  active  boundary-layer  control  system,  which  is  described  in  some 
detail  by  Saripalli  and  Simpson  (1980),  was  installed  on  the  non-test  walls 
of  the  test  section  to  inhibit  undesirable  flow  three-dimensionality  and  to 
prevent  separation.  Because  the  static  pressure  in  the  test  section  is  time 
varying  in  these  unsteady  experiments,  no  passive  boundary  layer  control  can 
be  used  that  depends  on  a  steady  test  section  pressure  higher  than  the 
pressure  outside  the  tunnel.  Highly  two-dimensional  wall  jets  of  high- 
velocity  air  were  introduced  at  the  beginning  of  each  of  the  three  sections 
of  the  tunnel.  At  the  later  two  streamwise  locations  the  incoming  boundary 
layer  was  partially  removed  by  a  highly  two-dimensional  suction  system. 

To  examine  the  two-dimensionality  of  the  mean  boundary  layer  flow, 
smoke  wa3  introduced  in  the  tunnel  and  a  sheet  of  laser  light  produced  by  a 
cylindrical  lens  was  used  to  illuminate  the  smoke.  Upstream  of  separation, 
negligible  spanwise  diffusion  of  the  smoke  was  observed,  indicating  no  gross 
flow  three-dimensionality.  Wall  jets  and  suction  boundary-layer  controls 
were  adjusted  to  produce  a  nearly  two-dimensional  flow  pattern  for  all  the 


three  flows. 
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2.4  Hot-Wire  Anemometer 

Miller  (1976)  type  integrated  circuit  hot-wire  anemometers  and 
linearizers,  as  modified  by  Simpson,  Heizer  and  Nasburg  (1979)  were  used  for 
the  present  experiments.  The  frequency  response  was  flat  up  to  7.5  kHz  for 
an  overheat  ratio  of  0.7.  This  moderately  high  overheat  ratio  was  used  for 
two  reasons.  First  as  shown  by  Wood  (1975)  the  range  of  flat  frequency 
response  is  improved  with  a  higher  overheat  ratio.  Secondly  to  minimize  the 
affect  of  air-temperature  oscillation  on  the  hot  wire  behavior.  A  standard 
TSI  model  1243-T.2  cross-wire  probe  was  used  for  measurements.  The  closest 
*:o  the  wall  that  this  probe  could  safely  make  measurements  was  about  0.035 
inch. 

The  traversing  mechanism  used  for  the  boundary  layer  velocity 
measurements  was  mounted  on  the  supporting  frame  for  the  upper  wall  except 
for  the  k  =  1.03  flow  with  the  "roof  damper".  In  that  case  the  hot  wire 
support  was  mounted  on  the  bottom  wall  and  the  hot-wire  stem  was  fed  through 
a  hole  in  the  test  wall  for  measurements  beyond  2.25  m.  A  cathetometer  was 
used  to  accurately  locate  the  probe  sensor  from  the  wall  with  an  uncertainty 
of  about  ±0.002  inches.  The  detailed  streamwise  free-stream  velocity 
distributions  were  obtained  using  a  model  TSI  121 8— T 1 . 5  model  probe  mounted 
on  a  mobile  cart  that  was  easily  positioned  along  the  flow.  Hot-wire 
calibrations  were  made  using  a  TSI  model  1125  calibrator.  A  standard  TSI 
model  1015C  correlator  was  used  to  obtain  difference  of  signals  from  two 
wires  of  the  cross-wire  probe. 

2.5  Laser  Anemometer 
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The  laser  anemometer  (LDA)  used  in  these  experiments  is  described  in 
some  detail  by  Simpson  and  Chew  (1979).  It  is  a  two-velocity-component 
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(U,V)  directionally  sensitive  fringe-type  system  that  has  been  by  Simpson  et. 
al.  (1983)  and  Simpson  and  Shivaprasad  (1983).  The  unshifted  and  21.505  MHz 
Bragg-cell  shifted  beams  lie  in  a  horizontal  plane  and  measure  the 
streamwise  velocity  with  vertical  fringes.  The  unshifted  and  -15.158  MHz 
Bragg  cell  shifted  beams  lie  in  a  vertical  plane  and  measure  (-V  cos  4.4°  + 
W  sin  4.4°)  with  almost  horizontal  fringes.  The  21.505  MHz  ana  -15.158 
MHz  beams  form  a  third  fringe  pattern  that  measures  (U  -  V  cos  4.4°  +  W  sin 
4.4°)  around  36.663  MHz.  Since 

u2  and  (-v  cos  4.4  +  w  sin  4.4)2 

were  measured  independently  and  uw  was  presumed  very  small,  the  Reynolds 
shearing  stress  -uv  resulted  from  this  measurement.  To  increase  the  signal- 

A 

to-noise  ratio  during  U  and  V  measurements,  vertical  and  horizontal  beams 
were  switched  off,  respectively,  and  the  central  beam  was  blocked  during  (U 

-  V  cos  4.4°  +  W  sin  4.4)//2  measurements.  Signal  processing  was  by  fast- 
sweep-rate  sampling  spectrum  analysis,  as  described  by  Simpson  and  Barr 
(1975). 

Difficulties  with  seeding  a  highly  turbulent  flow  are  well  known  since 
such  flows  are  characterized  by  intense  mixing.  This  problem  is  even  more 
aggravated  with  a  high  amplitude  unsteady  flow;  to  have  the  required  smoke 
concentration  at  each  phase  the  volumetric  flow  rate  needs  to  vary  with  the 
velocity.  The  smoke  is  generated  by  six  adjustable  Laskin  nozzles,  each  of 
which  blows  compressed  air  through  4  orifices  1  mm  in  diameter  into  dioctal 
phthalate  (DOP).  The  DOP  is  atomized  by  the  shearing  action  of  the 
compressed  air  jets  and  produces  a  mean  particle  diameter  of  about  1  pm. 
Here  0.057  cubic  meters  per  minute  of  smoke  at  a  concentration  of  about  4.8 

_  3 

x  10  kg  per  cubic  meter  of  blown  air  was  used.  The  smoke  was  distributed 
along  the  tunnel  width  just  upstream  of  the  test  section  by  33  evenly  spaced 
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holes  in  a  cylindrical  manifold,  which  is  supplied  with  the  smoke  into  both 
its  ends.  This  constant  volumetric  flow  rate  resulted  in  a  low  smoke 
concentration  and  consequently,  low  data  rate  during  the  accelerating 
velocity  phases  of  the  cycle  and  at  higher  velocity  phases  away  from  the 
wall,  where  velocities  are  higher. 

A  special  solenoid-controlled  butterfly  valve  which  provides  a 
variation  in  the  volumetric  flow  rate  was  designed  and  connected  between  the 
smoke  supply  and  uuo  end  of  the  smoke  manifold.  A  signal  from  a  pressure 
transducer  mounted  flush  with  the  tunnel  wall  in  the  upstream  wind  tunnel 
honeycomb  section  was  used  to  activate  the  solenoid.  An  electronic  circuit 
provided  a  time  delay  between  the  transducer  signal  and  activation  of  the 
solenoid  to  account  for  the  time  that  the  flow  takes  to  travel  from  the 
honeycomb  to  the  point  of  measurement.  The  valve  was  open  during 
accelerating  phases  of  the  cycle  and  closed  during  the  decelerating  phases. 
Trials  for  various  combinations  of  time  delays,  the  time  for  which  the  valve 
was  open,  and  the  smoke  flowrate  from  the  other  end  of  the  manifold  were  not 
very  successful  in  producing  a  more  uniform  data  rate  during  all  phases  of 
the  cycle.  However,  the  data  rate  was  improved  by  recirculating  the  used 
smoke  in  the  wind-tunnel  room  in  addition  to  that  obtained  from  the  smoke 
system. 

_  i  _  i 

The  minimum  data  rate  required  is  59S  and  95S  for  f  =  0.596  and 
0.954  Hz  flows  respectively.  But  since  new  signals  are  not  spaced  equally 
in  tim°,  a  higher  data  rate  is  necessary.  With  a  good  SNR  (above  15“20  dB) 
a  data  rate  of  400  new  signals  per  second  for  most  of  the  oscillation  cycles 


was  obtained. 
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2.6  Data  Reduction 

The  velocity  at  any  instant  can  be  written  as 

U(t)  =U+U+u=U+u 

where  0  is  the  time-averaged  velocity  U  is  the  periodic  velocity 
contribution  and  u  is  the  fluctuation  contribution. 

U  is  simple  time  averaged  of  the  velocity  value,  U  for  each  phase 
angle  is  obtained  by  ensemble-averaging  the  velocity  measurements  for  that 
phase  over  many  cycles.  Ensemble-averaged  mean  square  turbulence  quantities 
were  obtained  by  subtracting  the  smoothed  ensembled-averaged  velocity  for  a 
cycle  from  the  raw  data  and  averaging  the  square  over  many  cycles.  Mean  and 
turbulence  quantities  were  averaged  over  four  adjacent  bins  of  the  cycle  to 
produce  15°  averages  for  each  of  2^  bins. 

Results  were  obtained  by  averaging  over  200  oscillation  cycles. 
However,  for  laser  anemometer  measurements  at  some  phases  of  the  oscillation 
cycle  number  of  data  were  less  than  200,  as  data  carried  over  from  the 
previous  bins  were  discarded  during  data  processing.  Furthermore,  at  each 
phase  of  the  cycle  data  points  falling  beyond  three  standard  deviation  were 
discarded.  This  amounted  to  no  more  than  1-2?  rejection  of  data. 

Close  to  the  wall  and  up  to  distances  where  intermittent  backflow 
exists  at  any  phase  of  the  oscillation  cycle,  laser  anemometer  measurements 
were  made.  Cross-wire  anemometer  measurements  were  also  made  nearly  at  the 
same  axial  location.  To  increase  the  range  of  validity  of  X-wire  data,  in 
the  detached  flow  region,  the  X-wire  probe  was  inclined  10°  to  the  test 
wall,  i.e.,  mounted  parallel  to  the  top  wall  for  k  =  0.61  and  1.33  flows. 
In  data  reduction,  mean  and  turbulence  quantities  were  converted  back  into 
the  laboratory  coordinate  system. 
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Due  to  the  fact  that  the  laser  measurements  beyond  240  mm  from  wall 
were  not  possible  (restricted  by  movement  of  the  table  on  which  laser  and 
optics  were  mounted)  composite  profiles  of  mean  velocity  and  turbulence 
quantities  were  obtained  by  merging  the  data  files  obtained  by  these  two 
methods.  In  addition  to  providing  data  beyond  240  mm  from  the  test  wall, 
cross-wire  measurements  also  gave  some  region  of  overlap  with  LDA  data. 
Only  valid  X-wire  anemometer  data  were  retained  and  the  following  criteria 
were  used  to  check  the  validity: 

(a)  A  Gaussian  velocity  probability  distribution  is  approximately 
valid,  so  that  there  should  be  no  flow  reversal  (Simpson,  1976)  and  the  data 
are  valid  if 


y  ~  <  0.3  (1) 

U 

(b)  The  magnitude  of  the  flow  incidence  angle  to  the  cross-wire  bi¬ 
sector  is  less  than  or  equal  to  30°  then  the  data  are  considered  valid, 


i.e. , 
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By  the  above  criteria  the  cross-wire  data  cease  to  be  valid  as  the  backflow 
region  is  approached.  Valid  data  show  good  agreement  with  laser 
measurements  at  locations  where  both  measurements  are  available. 

Estimated  uncertainties  for  the  hot  wire  are  AU  =  ±3.2$,  AV2  and  Au2  = 
±10$,  A(-uv)  =  ±20$  and  for  the  laser  anemometer  are  AU  and  A V  =  ±0.25  m/s, 
Au2  and  Av2  =  ±4$  of  the  maximum  profile  value,  A(-uv)  =  ±20$  and  AX 

pu 

±0.02  for  all  phases.  For  the  phases  where  the  total  number  of  data  samples 
was  much  less  than  200  (after  discarding  the  sample  carried  over  from  the 
previous  bins)  the  uir  erta  inties  were  higher  and  the  LDA  data  were 
discarded . 
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3.  DESCRIPTION  OF  TEST  FLOWS 

The  ensembled-averaged  free-stream  velocity  Ug  outside  the  boundary 
layer  can  be  expressed  in  terms  of  its  Fourier  components  Ung  as 
„  ®  U 

U  =  0  (l  +  I  cos  (unt  -  <£_„)]}  (3) 

e  e  n=1  0 

e 

Figures  4  and  5  show  the  mean  free-stream  velocity  and  amplitude  of 
higher  harmonics  for  all  the  three  flows.  Upstream  of  2.8  m  the  first 


harmonic  ratio  - —  is  ~  0.8  for  k  =  0.61  and  1.33  flows  and  about  0.12  for 
°e 

k  =  1.03  flow  with  "roof  damper".  Second  harmonic  amplitude  to  mean  ratios 
are  8-10$,  15-20$  and  6-7$  whereas  third  harmonic  amplitude  to  mean  ratios 
are  4-6$,  3~7$  and  3~4$  of  the  mean  velocity,  respectively. 

Using  these  streamwise  distributions  the  free-stream  streamwise 
pressure  gradient  can  be  calculated  from  the  unsteady  Bernoulli  equation  and 
the  mean  and  first  harmonics  and  is  given  by 

1  -  -  [0  0  +  1  0.  0.  ]  +  [ (u  0;  +  0.  O')2 

p  dx  e  e  2  le  le  e  le  le  e 


+  Uf  (0  <t>’  -  U))2]  2  cos  (cot  -  <j>.  -  180°  -  Y.  ), 

le  e  le  le  le 


where 


y.  =  arctan 
le 


Ule(Vle  - 

0  u!  +0,0 

e  le  le  e 


Here  prime  denotes  streamwise  derivatives.  These  derivatives  are  evaluated 

at  a  given  streamwise  location  by  differentiating  the  least  square  curve  fit 

of  a  quadratic  model  of  the  five  data  points  nearest  to  that  location. 

Figures  6a-6c  show  the  non-dimensional  mean,  maximum  and  minimum 
dC 

pressure  gradients  along  the  centerline  of  the  test  wall.  Here  Cp  = 

2  ( P-P  )/pU2.  where  i  denotes  the  nominal  mean  free-stream  entrance 
l  ei 


condition.  Unlike  the  moderate  amplitude  flows  of  Simpson  et  al .  ( 1983a, b) 
in  the  diverging  section  a  small  favorable  pressure  gradient  exists  during 
the  high  velocity  phases  for  all  the  three  flows. 

Figures  7a-7c  show  the  phase  angle  of  the  first  harmonics  of  the  free- 

stream  velocity  <f>.  and  the  pressure  gradient  Y,  +  d>,  +  180  and  show  that 
le  le  le 

the  first  harmonic  pressure  gradient  strongly  lags  the  local  free-stream 
velocity  in  the  converging  section  of  the  tunnel  upstream  of  2  m.  The  lag 
is  considerably  lower  in  the  diverging  section.  After  detachment  the 
oscillating  pressure  gradient  only  slightly  lags  (or  leads)  the  velocity 
oscillation  with  the  onset  of  pressure  gradient  relaxation,  as  also  observed 
by  Simpson  et  al .  (1983a, b)  in  moderate  amplitude  unsteady  flows. 

Free-stream  velocity  distribution  along  the  tunnel  length  (Figures  8- 
10)  show  a  disturbance  travelling  downstream.  The  changes  in  free  stream 
velocities  are  confirmed  by  X-wire  and  normal  wire  measurements.  Although 
the  effect  with  high  frequency  flows  (k  =  1.33  and  k  =  1.03  with  "roof 
damper")  is  quite  pronounced,  a  weaker  affect  is  evident  with  k  =  0.61  flow. 
This  disturbance  has  a  phase  shift  of  about  ir/3  for  X  =  0.61  m  to  X  =  4.27  m 
for  k  =  1.33  flow.  For  k  =  1.03  flow  with  "roof  damper",  this  effect  is 
seen  only  beyond  3  m  with  a  phase  shift  of  about  tt / 4  from  X  =  3.30  m  to  X  = 
4.37  m.  This  effect  was  not  observed  in  moderate  amplitude  flows  of  Simpson 
et  al . 

Study  of  the  velocity  and  pressure  gradient  oscillations  shows  that 
upstream  of  the  throat  where  the  maximum  velocity  increases  as  the  tunnel 
area  decreases  there  are  some  phases  where  an  adverse  pressure  gradient 
exists  as  shown  in  Figures  11-13  for  all  the  three  flows.  In  the  diverging 
section  of  the  tunnel,  where  large  backflow  exists,  an  adverse  pressure 
gradient  is  present  for  most  of  the  phases  of  the  cycle  and  higher  harmonics 
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of  the  pressure  gradient  are  produced  that  subject  the  boundary  layer  to 
several  pressure  gradient  oscillations  per  cycle  of  the  velocity.  For  k  = 
1.33  flow,  higher  harmonics  of  the  pressure  gradient  are  produced  upstream 
of  detachment  also. 


4.  EXPERIMENTAL  RESULTS 

Ensemble-averaged  velocity  profiles  were  obtained  using  cross-wire  and 
laser  anemometers.  In  regions  where  both  the  laser  and  cross-wire 
anemometer  produce  valid  data,  the  results  agreed  within  the  experimental 
uncertainties. 

Ensemble-averaged  velocity  data  for  various  phases  at  X  =  0.53  m  and 
1.45  m  for  large  amplitude  flows  and  at  X  =  1.63  m  for  k  =  1.03  flow  with 
"roof  damper"  are  shown  in  Figures  14  to  16.  For  the  k  =  1.03  flow  with 
"roof  damper"  at  all  the  phases  and  at  the  high  free-stream  velocity  phases, 
i.e.,  45°  <  uit  S  3 1 5 0  ,  of  the  k  =  0.61  and  1.33  flow  the  shape  of  the 
velocity  profiles  are  essentially  the  same  as  observed  by  Simpson  et  al . 
( 1983a, b)  for  the  steady  and  moderate  amplitude  flow  and  are  similar  to  that 
for  low  Reynolds  number  attached  turbulent  boundary  layers.  A  semi- 
logarithmic  region  is  clearly  evident. 


/u2/Ue  and  /v2/U  collapse  on  top  of  another  for  these  high  velocity 
and  accelerating  phases,  as  shown  in  Figures  17-18.  Figures  19a  and  19b 
give  law-of-the-wall  plots  for  these  phases.  The  mean  wall  shear  stress 
values  were  obtained  by  adopting  Coles  and  Hirst  (1969)  method  which 


requires  that  for  a  given  U  ,  U  =  16.23  at  y 
of-the-wall  equation 


A  +  1  + 

U  =  -  in  y  +5.0 


100.  The  logarithmic  law- 


/s'1  v.v  v;.>  v,  .v.-.  a. 


(5) 


where  <  is  0.41  ,  is  also  shown  in  Figures  19a- 19b  by  a  solid  line.  Skin 
friction  coefficient  values  obtained  by  this  method  are  presented  in  tables 


C,  -  0.2-6  x  ,0-°-678H  x  Se--°-268 
f  8 


1~3  (Appendix  A).  As  can  be  expected  for  these  flows  with  adverse  pressure 
gradients,  the  skin  friction  decreases  in  the  streamwise  direction  until  the 
reattachment  process  begins.  The  skin  friction  coefficient  was  also 
calculated  using  the  Ludwieg-Tillmann  (1950)  relation  (which  is  valid  for 
the  phases  with  the  attached  flow), 

Q 

(6) 

The  skin  friction  coefficients  obtained  by  these  two  methods  show  good 

agreement,  within  the  experimental  uncertainties,  as  shown  in  Figure  20. 

Variation  of  6*,  0,  H  and  Ret  for  all  the  three  flows  before  detachment 

0 

are  shown  in  Figures  21a-21d,  22a-22d  and  23a. 

At  low  velocity  phases  for  the  k  =  0.61  and  1.33  flows  close  to  the 
wall  intermittent  backflow  exists  upstream  of  the  throat  as  shown  in  Figures 
24-25  with  the  fraction  of  time  that  there  is  downstream  flow  Ypu  less  than 
unity.  These  values  were  obtained  from  direct  LDA  measurements  and  also 

A  A 

from  the  LDA,  U  and  u2  data  using  the  Gaussian  probability  distribution 
(Simpson  1976) 
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In  two  moderate  amplitude  flows  studied  by  Simpson  et  al .  (  1  983a,  b)  this 
intermittent  backflow  upstream  of  detachment  was  not  observed.  Ensemble- 
averaged  velocity  profiles  at  these  phases,  i.e.,  315°  £  u>t  S  45°  (as  shown 
in  Figures  14  and  15  for  15°,  30°,  345°  and  360°)  are  quite  different. 
These  velocity  profiles  at  the  low  velocity  phases  where  the  flow  is 
accelerating,  i.e.,  wt  S  45°  are  very  similar  to  the  strongly-accelerated 


sink  flow  studied  by  Simpson  and  Wallace  (1975).  The  unsteady  flow 
acceleration  parameter 
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has  large  values  at  these  phases  as  shown  in  Figures  27“28  (but  not  for  k  = 
1.03  flow  as  shown  in  Figure  29). 

Simpson  and  Wallace  (1975)  observed  that  for  a  steady  flow  with  K  >  3.6 
_  6 

x  10  at  low  momentum  Reynolds  numbers,  relaminar ization  of  an  initially 

~  _  6 

turbulent  boundary  layer  could  occur.  For  the  k  =  0.61  flow,  K  >  3.6  x  10 

over  20°  of  the  cycle  during  rapid  acceleration  from  the  minimum  velocity 

phases.  For  the  same  phases  in  the  k  =  1.33  flow,  K  is  over  10  times 

A  A 

greater.  In  both  these  flows  nearly  zero  correlation  coefficients  -uv/uv 
were  observed  at  these  phases  (Figure  30).  These  large  K  values  are 
sufficient  to  produce  laminarescent  ensemble-averaged  profiles  during  a 
short  portion  of  the  cycle,  but  are  not  imposed  long  enough  to  produce  true 
relaminarization.  Velocity  profiles  at  the  locations  (Figures  14  and  15) 
for  (ut  =  3^5°  and  360°,  when  the  flow  is  highly  decelerating  and  similar  to 
typically  high  adverse  pressure  gradient  flows  (Simpson  1981  a). 

Upstream  of  detachment  the  ensemble-averaged  velocity  profiles  exhibit 
a  negligible  phase  angle  variation  through  the  boundary  layer  for  all  three 
flows  (Figures  3 1  ~  3  3 )  -  Simpson  et  al .  (  1983b)  showed  that  the  phase  angle 
variation  of  the  ensemble-averaged  oscillation  is  zero  within  a  semi- 
logarithmic  region.  The  ensemble-averaged  turbulence  profiles  for  the 
phases  without  flow  reversal  are  similar  to  those  for  the  comparable  steady 
flow  and  u2  and  v2  are  in  phase  with  U.  Near  detachment,  outer  region 
results  for  the  ensemble-averaged  velocity  profiles  at  the  phases  without 
flow  reversal  follow  the  universal  velocity  defect  law  proposed  by  Schofield 
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(1980).  Ensemble-averaged  velocity  profiles  in  the  detaching  and 
reattaching  parts  of  the  flows  are  presented  in  Figures  34-43.  Velocity 

U  V 

profiles  for  the  "detaching"  phases  are  presented  as  s—  vs  x  ,  whereas  for 

Ue  6 

the  reattaching  phases  absolute  velocities  are  presented  (this  is  because  of 
absence  of  the  core  for  these  low  velocity  phases).  For  k  =  1.03  flow  with 
"roof  damper"  all  the  velocity  profiles  (except  at  X  =  4.12  m)  are  presented 
U  y 

as  x  vs  x  .  Also  for  this  flow,  as  evident  from  Figures  40-43,  velocity 

U  6 

e 

profiles  almost  collapse  over  one  another  except  at  the  phases  where  the 
velocity  profile  shows  rapid  changes  because  of  the  travelling  wave. 

For  the  k  =  0.61  and  1.33  flows,  backflow  velocities  are  very  large  and 
for  the  k  =  1.33  flow  at  some  phases  are  even  larger  in  magnitude  than  the 
free-stream  velocity.  Figures  34,  37  and  40  show  distinct  semi-logarithmic 
velocity  profiles  for  the  part  of  the  cycle  as  soon  as  there  is  no  mean 
backflow.  Figures  44  and  45  give  law-of-the-wall  plots  for  these  phases  for 
the  k  =  0.61  and  1.33  flows.  Values  of  wall  shear  stress  deduced  from  the 
attached  flow  logarithmic  "law  of  the  wall"  velocity  profiles  are  in  good 
agreement  with  obtained  using  Ludwieg  and  Tillmann  equation  as  shown  in 
Figures  46  and  are  also  within  experimental  uncertainty  of  the  measured 
Reynolds  shear  stresses  in  the  semi-log  region.  This  reattachment  of  the 

3  U 

flow  is  apparently  due  to  the  increasing  rate  of  strain  ,  the  outer 
region  shear  stress  for  the  thinning  shear  layer,  and  the  convection  of  the 
upstream  flow.  Note  that  the  semi-log  region  velocities  remain  almost 
independent  of  time  for  105°  <  uit  <  165°. 

During  larger  backflows,  significant  turbulence  production  occurs  near 
the  wall  and  all  the  Reynolds-averaged  momentum  equation  terms 


«  V  V.V. -W.V„  J 


31311  C9) 

at  8x  Sy  p  9x  ay  1,1 

are  important.  Figures  47-50  show  that  decreasing  Reynolds  shearing  stress 

and  correlation  coefficients  of  -0.5  or  so  are  produced  near  the  wall  during 

these  larger  backflow  phases.  Typically  large  correlation  coefficients  of 

0.5  are  in  the  outer  region  where  most  of  the  turbulence  energy  production 

and  the  maximum  shearing  stress  occur.  During  the  reattaching  or  "washout" 

of  detached  flow  phases,  the  backflow  velocities  decrease  and  the  near  wall 

shearing  stresses  and  correlation  coefficients  increase  toward  attached 

u 2  v  2 

flow  values.  Figures  51-56  show  ensemble-averaged  jrz  and  jj?  profiles  at 

e  e 

some  locations  in  detached  flow. 

As  suggested  by  Simpson  et  al .  (1981)  the  near  wall  backflow  is 

dominated  by  viscous  affects  because  the  Reynolds  shearing  stress  is 

relatively  unimportant  and  the  mean  flow  can  be  considered  to  be  uncoupled 

U 

for  the  turbulence.  A  velocity  profile  model  of  the  type  p-  =  f(y/N)  fits 

N 

the  flow  structure.  For  the  k  =  0.61  and  1.33  flows  a  semi -logar  i  thmi  c 
overlap  region  described  by 

-J —  -  -0.3  tn|i|  *  B  (10) 

M  " 

A 

seem  to  exist  for  0.02  <  y/N  <  0.15,  where  B  is  a  constant.  For  the  k  = 
1.03  flow  with  "roof  damper"  this  semi-logarithmic  region  also  exist  (as 
shown  in  Figure  57)  but  backflow  thicknesses  are  very  large  (>25  cm)  and  it 
appears  there  is  some  additional  length  scale  which  influences  the  outer 
part  of  the  backflow.  We  are  continuing  to  investigate  this  behavior. 

Ensemble-averaged  velocities  in  the  backflow  lead  the  free-stream 
velocity  by  a  large  amount.  Similar  to  the  moderate  amplitude  flows, 
(Simpson  et  al .  ,  1983),  the  phase  angle  of  first  harmonic  does  not  vary 
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appreciably  with  y  in  the  nearest  wall  region.  Far  downstream  where  the 
ensemble-averaged  velocity  is  always  negative,  <j>  1  ,  and  are  nearly 

independent  of  y  for  y  <  N.  Thus  it  appears  that  the  ensemble-averaged 
flow  near  the  wall  behaves  like  a  quasi-steady  flow  when  normalized  on 
and  N. 

Figures  31 “33  show  that  the  ensemble-averaged  velocity  in  the  backflow 
has  a  progressively  greater  phase  lead  as  one  proceeds  downstream.  This 
phase  lead  with  respect  to  the  local  freestream  is  higher  for  the  k  =  0.61 
flow  compared  to  the  k  =  1.33  flow  and  is  even  higher  for  the  k  =  1.03  flow 

with  "roof  damper".  Upstream  of  detachment  the  Reynolds  shearing  stresses 

~  ~  ~ 

-uv  and  the  normal  stresses  u  and  v2  are  in  phase  with  U.  Downstream  of 
detachment  u2  and  v2  are  substantially  out  of  phase  with  -uv  and  -uv 
slightly  lags  Uj  (about  20°)  as  shown  in  Figures  58-60.  Hysteresis  in 
relationship  in  mean  velocity  and  turbulence  parameter  have  been  observed  as 
in  the  case  of  moderate  amplitude  flows  studied  by  Simpson  et  al .  Figures 
61-63  show  IU..I/U  vs.  Y  .  for  all  the  three  flows.  The  hysteresis 
loops  are  clockwise  (as  in  moderate  amplitude  flows)  and  are  progressively 
greater  in  the  downstream  direction  (for  the  larger  amplitude  flows,  loops 
are  not  complete  because  of  absence  of  core  at  some  phases).  Data  are 
scattered  at  lower  values  of  Ypu,  because  of  the  large  percentage 
uncertainty  when  Y  is  near  zero.  Hysteresis  loops  are  larger  for  the  k  = 

A  A 

1.33  flow.  In  all  the  plots  |U^|/Ue  increases  substantially  during  the 


minimum  Y  .  phases, 
pu  min  K 


In  Figures  64  and  66  data  are  presented  for 
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vs.  -*•  for  the  k  =  1.03  flow  with  "roof  damper".  Data  for  large  amplitude 
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flow  (k  =  1.33)  are  presented  on  another  parameter  cross-plot  i.e., 


in  Figure  65.  Similar  hysteresis  loops  are  observed  except  for 


/u2  /v2 


the  k  =  1.03  flow  with  "roof  damper"  where  double  loops  exist.  We  are 
continuing  to  investigate  this  phenomenon.  For  the  k  =  1.03  flow  with  "roof 
damper"  at  some  phase  there  appears  to  be  no  hysteresis  because  of  a  nearly 


constant  value  of  the  shape  parameter,  H. 

Figure  67  presents  the  normalized  maximum  backflow  velocity  |U„ I /U  vs. 

1  N 1  e 

1  for  all  the  three  flows.  The  data  for  large  amplitude  flows  fall  between 
H 

the  steady  separating  flow  line  (solid  line)  of  Simpson  and  Shivaprasad 


(1983), 


M, 


0.807  (1  -  -)  -  0.577 
H 


and  that  given  by  (dashed-line) 
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which  was  derived  (Simpson  and  Shivaprasad,  1983)  from  an  approximate 

velocity  profile  based  upon  the  assumption  that  ensemble-averaged  profiles 

for  phases  when  97  .  / 3 1  £  0  are  rounded  near  the  wall  and  near  the  free- 

pu  mm 

stream  for  a  given  |UN|/Ug.  As  in  the  moderate  amplitude  flows,  |U^|/Ue 

reaches  large  values  with  negative  values  of  x  .  As  shown  by  Simpson  and 

H 

Shivaprasad  (1983),  for  phases  of  a  cycle  when  9Ypu/9t  is  positive 
different  hysteresis  loops  occurs  on  these  plots,  indicating  a  different 
velocity  profile  and  turbulence  structure  from  increasingly  separated 


5.  SUMMARY  OF  RESULTS 


It  is  clear  from  these  measurements  and  the  previous  studies  of  Simpson 
et  al.  (1983)  and  Simpson  and  Shivaprasad  (1983)  in  moderate  amplitude 
flows,  that  the  oscillation  waveform  and  amplitude  strongly  influence  the 
detached  flow  behavior.  This  is  an  important  result  since  practical 
unsteady  flows  have  non-sinusoidal  waveforms  and  significantly  large 
amplitudes.  Prior  to  the  current  work  no  detailed  velocity  and  turbulence 
measurements  on  large  amplitude  type  flows  had  been  made  even  though 
practical  unsteady  flows  in  turbomachines  have  non-sinusoidal  waveforms  and 
significantly  large  amplitudes. 

During  the  phases  that  the  flow  is  accelerating,  intermittent  backflow 
occurs  only  near  the  wall  in  the  rapidly  diverging  portion  of  the  test 
channel.  For  large  amplitude  flows  during  decelerating  and  minimum  velocity 
phases,  backflow  occurs  in  the  diffuser  section  and  near  the  wall  in  the 
converging  section  as  well.  This  behavior  was  not  observed  by  Simpson  et 
al .  in  moderate  amplitude  flows.  This  is  an  interesting  result  since  it 
could  be  used  as  a  possible  means  of  increasing  near  wall  mixing  of 
turbulent  fluid  to  enhance  heat  and  mass  transfer. 

For  all  the  flows  studied  higher  harmonics  of  the  pressure  gradient  are 
produced  that  subject  to  the  boundary  layer  to  many  pressure  gradient 
oscillations  for  cycle  of  the  velocity. 

Upstream  of  detachment  the  ensemble-averaged  velocity  profiles  exhibit 
a  negligible  phase  angle  variation  through  the  boundary  layer  and  u2,  v2  and 
-uv  are  in  phase  with  U.  The  velocity  and  turbulence  in  the  absence  of  flow 
reversal  behaves  in  a  manner  as  the  steady  flow  and  the  flow  is  quasi- 
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Downstream  of  detachment  large  amplitude  and  phase  variations  develop 
through  the  flow  and  there  are  non-quasi-steady  affects  on  the  ensemble- 
averaged  flow  structure.  Backflow  velocities  even  larger  than  the  free- 
stream  have  been  observed.  Near  the  wall  in  the  backflow  region,  the 
ensemble-averaged  velocity  leads  the  free-stream  velocity  by  a  large  amount. 
The  turbulence  structure  progressively  lags  the  ensemble-averaged  flow 
oscillations  with  -uv  lagging  Uj  in  the  backflow  by  about  20°,  as  in  the 
moderate  amplitude  flows.  The  ensemble-averaged  backflow  near  the  wall 
behaves  like  a  quasi-steady  flow  when  normalized  on  and  N.  Study  of  the 
ensemble-averaged  backflow  velocities  for  the  k  =  1.03  flow  with  "roof 
damper"  and  very  large  backflow  thicknesses  suggest  that  there  are  some 
other  length  scales  for  the  outer  part  of  the  backflow  region.  We  are 
continuing  to  investigate  this  behavior. 

Low  and  negative  values  of  Reynolds  shearing  stress  and  correlation 
coefficient  during  phases  with  flow  reversal  have  been  measured.  This 
indicates  that  much  of  the  reversed  flow  turbulence  is  left  over  from 
earlier  phases  when  the  turbulence  production  was  large  and  that  the  adverse 
pressure  gradients  are  large  enough  tc  produce  significant  backflow  with 
negative  shearing  stress.  This  is  a  non-quasi-steady  effect  which  the 
current  local  field  turbulence  models  can  not  account  for.  In  the  detached 


and  reattaching  flow,  ensemble-averaged  velocity  profiles  and  turbulence 
structure  show  the  hysteresis  that  occur  in  detaching  and  reattaching  flows, 
which  is  again  a  non-quasi-steady  behavior. 


Additonal  discussion  of  these  data  and  their  meaning  will  be  contained 
in  Agarwal  and  Simpson  ( 1988a, b). 
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6.  PERSONNEL 

During  1984-1986,  Mr.  S.  Ahn  made  single  hot-wire  anemometer 
measurements  for  the  k  =  0.61  flow,  resulting  in  his  M.S.  degree  thesis 
(Ahn,  1986).  Dr.  Y.-T.  Chew,  who  is  senior  lecturer  at  the  National 
University  of  Singapore,  worked  on  multiple  hot-wire  anemometer  measurement 
techniques  during  his  July  1984  -  February  1985  visit  to  VPI&SU,  resulting 
in  two  papers  (Chew  and  Simpson,  1988a, b).  Dr.  N.K.  Agarwal  joined  this 
research  group  in  May  1985  after  completing  his  Ph.D.  work  at  the  University 
of  Adelaide,  Australia.  In  addition  to  this  report,  he  is  coauthor  of 
several  papers  on  this  research:  Agarwal  and  Simpson  (1986,  1987a,  1987b, 
1988a,  1988b). 

7.  INTERACTIONS 

A  number  of  oral  presentations  of  these  research  results  have  been 
made,  in  most  cases  accompanying  a  written  paper.  Simpson  (1985)  and 
Agarwal  and  Simpson  (1986,  1987a,  1987b)  were  presented  at  the  respective 
conferences  indicated  in  the  references.  Agarwal  and  Simpson  (1987a)  was 
also  presented  at  the  1986  American  Physical  Society  Fluid  Dynamics  Meeting 
in  Columbus,  Ohio. 
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Figure  6a.  Mean  pressure  gradient  distribution  along  the  tunnel  f- 

centerline  and  Jimits  of  the  unsteady  pressure  gradient,  k  =  0.61  flow; 

Cp  =  2 (P  —  P)!pUl  .  (  +  .maximum;  □  .mean;  *  .minimum). 
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Figure  6b.  Mean  pressure  gradient  distribution  along  the  tunnel 
centerline  and  limits  of  the  unsteady  pressure  gradient.  k  =  1.33  flow; 
Cp  =  2(P  —  P,)/ pUl,  .  (  •*•  .maximum;  -f  .mean;  q  .minimum). 


Figure  6c.  Mean  pressure  gradient  distribution  along  the  tunnel 
centerline  and  limits  of  the  unsteady  pressure  gradient,  k=1.03  flow 
with  'roof  damper';  Cp  =  2 (P  -  P,)lpUl  ;  (  a  .maximum;  x  .mean;  □ 
.minimum). 
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Figure  8a.  Ensemble-averaged  free-stream  velocity  along  the  tunnel 
length,  k  =  0.61  flow;  0.09  m  <  X  <  1.46  m  ; 

(  X  =  0.09m;  +  =  0.24m;  O  =  0.39m  ;  *  =  0.55m;  El  —  0.70m  ; 
A  -  0.85m;  ^  =  1.00m  ;  O  =  1.16m;  □=  1.31m  ;  ■=  1.46m). 
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Figure  8b.  Ensemble-averaged  free-stream  velocity  along  the  tunnel 
length,  k  =  0.61  flow.  1.61  m  ^  X  <  4.66  m  ;  (  A  =  1.61m; 
A  =  1.76m;  O  =  1.92m  ;  ♦  =  2.07m;  v  =  2.22m  ;▼  =  2.37m; 
O  =  2.53m  ;•  =  2.68m;  Cf  =  2.83m  ;  ■  =  2.98m;  ;□  =  3.14m.  ; 
W  =  3.29m;  XX  =  3.44 m  =  3.59m;  A  =  3.75m  ;  A  =  3.90m; 
<$  =  4.05m  ;  ♦  =  4.20m;  y  =  4.35m  ;  y  =  4.51m;  <=>  =  4.66m  ) 
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Figure  9a.  Ensemble-averaged  free-stream  velocity  along  the  tunnel 
length,  k  =  1 .33  flow.  0.15  m  <  X  <  1,52  m  ;  (X  =  0.15m; 
+  =  0.30nr;  O  -  0.46m  ;  *  —  0.61m;  [3  —  0.76m  ;  A  =  0.91m; 
v=  1.07m;  G=  1.22m;  □=  1.37m;  ■=  1.52m). 
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Figure  10b.  Ensemble-averaged  free-stream  velocity  along  the  tunnel 
lei  -th,  k  =  1.03  flow  with  'roof  damper'.  1.61m  <  X  <  4.66m  :  Legend 
same  as  in  figure  8b. 
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Figure  11a.  Free-stream  velocity  and  pressure  gradient  at  X  =  0.55  m, 
k  =  0.61  flow.  (  □  ,free  -  stream ;  a  , pressure  gradient ) 
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Figure  11b.  Free-stream  velocity  and  pressure  gradient  at  X=1.46  rn, 
k  =  0.61  flow.;  Legend  same  as  in  figure  11a, 
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Figure  11c.  Free-stream  velocity  and  pressure  gradient  at  X  =  2.22  m, 
k  =  0.61  flow.;  Legend  same  as  in  figure  11a. 
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Figure  lid.  Free-stream  velocity  and  pressure  gradient  at  X  =  2.83  m, 
k  =  0.61  flow;  Legend  same  as  in  figure  11a. 
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Figure  Ilf.  Free-stream  velocity  and  pressure  gradient  at  X  =  4.20 
k  =  0.61  flow;  Legend  same  as  in  figure  11a. 
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Figure  13a.  Free-stream  velocity  and  pressure  gradient  at  X=1.61  m, 
k=  1.03  flow  with  'roof  damper';  Legend  same  as  in  figure  11a. 
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Figure  13b.  Free-stream  velocity  and  pressure  gradient  at  X  =  2.98  ra, 
k=  1.03  flow  with  'roof  damper';  Legend  same  as  in  figure  11a. 
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Figure  13c.  Free-stream  velocity  and  pressure  gradient  at  X  =  3.29 
k=  103  flow  with  roof  damper';  Legend  same  as  in  figure  11a. 
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Figure  13d.  Free-stream  velocity  and  pressure  gradient  at  X  =  3.59  m 
k=  1.03  flow  with  'roof  damper';  Legend  same  as  in  figure  1  la. 
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Figure  13e.  Free-stream  velocity  and  pressure  gradient  at  X  =  4.51 
K=  1 .03  flow  with  'roof  damper';  Legend  same  as  in  figure  11a. 
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50°;  v  =  165°  ;  ■  =  180°)  open  symbols,  laser  anemometer; 
mbols,  cross-wire  measurements.  Solid  line  given  by  equation 


Figure  20.  Correlation  between  C,  obtained  from  Semi-Log  Region  and 
C,  obtained  from  Ludwieg-Tillman  equation  6;  90°  ^  cot  <,  270°; 
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Figure  21a.  Variation  of  displacement  thickness,  momentum  thickness, 
shape  factor  and  Reynolds  number  with  phase  angle  at  X  =  0.53  m, 

k  =  0.61  flow.  (  0  =  6‘,  X  =  6,  A  =  H,  +  =  RJ 
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Figure  21b.  Variation  of  displacement  thickness,  momentum  thickness, 
shape  factor  and  Reynolds  number  with  phase  angle  at  X=1.45  m, 
k  =  0.61  flow.  Legend  same  as  in  figure  21a. 
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Figure  21c.  Variation  of  displacement  thickness,  momentum  thickness 
shape  factor  and  Reynolds  number  with  phase  angle  at  X=1.83  m 
k  =  0.61  flow.  Legend  same  as  in  figure  21a. 
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Figure  21  d.  Variation  of  displacement  thickness,  momentum  thickness, 
shape  factor  and  Reynolds  number  with  phase  angle  at  X  =  2.20  m, 
k  =  0.61  flow.  Legend  same  as  in  figure  21a. 
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Figure  21e.  Variation  of  displacement  thickness,  momentum  thickness 
shape  factor  and  Reynolds  number  with  phase  angle  at  X  =  2.95  m 
k  =  0.61  flow.  Legend  same  as  in  figure  21a. 
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Figure  2 1  f .  Variation  of  displacement  thickness,  momentum  thickness, 
shape  factor  and  Reynolds  number  with  phase  angle  at  X  =  3.51  m, 
k  =  0.61  flow.  Legend  same  as  in  f'gure  21a. 
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Figure  21g.  Variation  of  displacement  thickness,  momentum  thickness, 
shape  factor  and  Reynolds  number  with  phase  angle  at  X  =  4.27  m, 
k  =  0.61  flow.  Legend  same  as  in  figure  21a. 
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Figure  22a.  Variation  of  displacement  thickness,  momentum  thickness, 
shape  factor  and  Reynolds  number  with  phase  angle  at  X  =  0.53  m, 
k=  1.33  flow.  Legend  same  as  in  figure  21a. 
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Figure  22b.  Variation  of  displacement  thickness,  momentum  thickness 
shape  factor  and  Reynolds  number  with  phase  angle  at  X^l.45  m 
k  =  1 .33  flow.  Legend  same  as  in  figure  21a. 
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Figure  22c.  Variation  of  displacement  thickness,  momentum  thickness, 
shape  factor  and  Reynolds  number  with  phase  angle  at  X=1.83  m, 
k=  1.33  flow.  Legend  same  as  in  figure  21a. 
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Figure  22d.  Variation  of  displacement  thickness,  momentum  thickness, 
shape  factor  and  Reynolds  number  with  phase  angle  at  X=:2.20  m, 
k=  1 .33  flow.  Legend  same  as  in  figure  21a. 
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Figure  22e.  Variation  of  displacement  thickness,  momentum  thickness, 
shape  factor  and  Reynolds  number  with  phase  angle  at  X  =  2.95  m, 
k=  1.33  flow.  Legend  same  as  in  figure  21a. 
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Figure  22f.  Variation  of  displacement  thickness,  momentum  thickness, 
shape  factor  and  Reynolds  number  with  phase  angle  at  X  =  3.51  m, 
k=  1 .33  flow.  Legend  same  as  in  figure  21a. 
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Figure  22g.  Variation  of  displacement  thickness,  momentum  thickness, 
shape  factor  and  Reynolds  number  with  phase  angle  at  X  =  4.27  m, 
k  =  1.33  flow.  Legend  same  as  in  figure  21a. 
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Figure  23a.  Variation  of  displacement  thickness,  momentum  thickness, 
shape  factor  and  Reynolds  number  with  phase  angle  at  X=1.63  m, 
k=  1.03  flow  with  'roof  damper'.  Legend  same  as  in  figure  21a. 
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Figure  23b.  Variation  of  displacement  thickness,  momentum  thickness, 
shape  factor  and  Reynolds  number  with  phase  angle  at  X  =  3.01  m, 
k=  1.03  flow  with  'roof  damper'.  Legend  same  as  in  figure  21a. 
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Figure  23c.  Variation  of  displacement  thickness,  momentum  thickness, 
shape  factor  and  Reynolds  number  with  phase  angle  at  X  =  3.28  m, 
k=  1.03  flow  with  'roof  damper'.  Legend  same  as  in  figure  21a. 
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Figure  23d.  Variation  of  displacement  thickness,  momentum  thickness, 
shape  factor  and  Reynolds  number  with  phase  angle  at  X  =  3.58  m, 
k=  1 .03  flow  with  'roof  damper'.  Legend  same  as  in  figure  21a. 


Figure  24b.  Variation  of  you  vs.  X  for 
different  phases  of  a  cycle,  k  =  0.61 
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Figure  24a.  Variation  of  ypu  vs.  X  for 
different  phases  of  a  cycle,  k  =  0.61 

flow;  ( - =  15°  ; - =  30°; 
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Figure  24c.  Variation  of  v.„  vs.  X  for 
different  phases  of  a  cycle.  k  =  0  61 
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Figure  24d.  Variation  of  .  vs.  X  for 
different  phases  of  a  cycle,  k  =  0.6 1 

•  low;  ! -  235  . - -  30C 

--dip  •  -  ~  330  . 

—  345 


360  l. 


Figure  25a.  Variation  of  ypu  vs.  X  for 
different  phases  of  a  cycle,  K ~  1 .33 
flow;  Legend  same  as  in  figure  24a. 


Figure  25b.  Variation  of  ypu  vs.  X  for 
different  phases  of  a  cycle,  k  =  1.33 
flow;  Legend  same  as  in  figure  24b. 
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Figure  25c.  Variation  of  vs.  X  for 
different  phases  of  a  cycle,  k  =  1 .33 
flow;  Legend  same  as  in  figure  24c. 


Figure  25d.  Variation  of  vs.  X  for 
different  phases  ••  a  cycle,  k  =  1 .33 
flow;  Legend  same  as  in  figure  24d. 


Figure  26a.  Variation  of  ypu  vs.  X  for 
different  phases  of  a  cycle,  k  =  1 .03  flow 
with  'roof  damper7;  Legend  same  as  in 
figure  24a. 


Figure  26b.  Variation  of  yBU  vs.  X  for 
different  phases  of  a  cycle,  k  =  1 ,03  flow 
with  'roof  damper';  Legend  same  as  in 
figure  24b. 
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Figure  26c.  Variation  of  ypu  vs.  X  for 
different  phases  of  a  cycle,  k  =  1 .03  flow 
with  roof  damper';  Legend  same  as  in 
figure  24c. 


Figure  26d.  Variation  of  vs.  X  for 
different  phases  of  a  cycle.  k=  1.03  flow 
with  roof  damper';  Legend  same  as  in 
figure  24d. 
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Figure  28 a.  Variation  of  acceleration  parameter,  K  with  phase  angle  at 
X  =  0.53  m,  k=  1.33  flow. 


Variation  of  phase  angle  of  first  harmonic 
nd  same  as  in  figure  31. 
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Figure  42a.  Ensemble-averaged  velocity  profiles,  U/U „  vs.  log)0(y/4),  at 
X  =  3.59  m,  k  =  1.03  flow  with  'roof  damper'  at  15°  <,  wt  <,  180°;  Legend 
same  as  in  figure  40a. 
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Figure  47a.  Ensemble-averaged  Reynolds  shearing  stress  —  uv/U l 
profiles  at  X  =  4.27  m,  k  =  0.61  flow  at  phases  120  °<  a>f^195°;  ( 
□  =  120°;  A  =  135°;  O  =  150°  ;  v  =  165°;  O  =  180°  ; 
Lf  =  195°  );  open  symbols,  laser  anemometer;  solid  symbols, 
cross-wire  measurements. 


Figure  47b.  Ensemble-averaged  Reynolds  shearing  stress  -  uv/U* 
profiles  at  X  =  4.27  m,  k  =  0.61  flow  at  phases  210  cot  <,285°:  ( 
□  =  210”;  A  =  225°;  o  =  240°  ;  v  =  255°;  O  =  270°  ; 
Cf  =  285°);  open  symbols,  laser  anemometer;  solid  symbols, 
cross-wire  measurements. 
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Ensemble-averaged  shear  stress  correlation  profiles  at  X 
=  1.33  flow  at  phases  195  °  <,  cot  <L  330°;  Legend  same  as  in 
ines  for  visual  aid  only. 
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Figure  50c.  Ensemble-averaged  Reynolds  shearing  stress  -  uvjUl 
profiles  at  X  =  3.59  m,  k  =  1.03  flow  with  'roof  damper'  at  phases  195 
°  <  cot  <,  270°;  Legend  same  as  in  figure  40a. 
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Figure  50d.  Ensemble-averaged  Reynolds  shearing  stress  —  uvjUl 
profiles  at  X  =  3.59  m,  k=1.03  flow  with  'roof  damper7  at  phases  285 
°  cot  <,  360°;  Legend  same  as  in  figure  40a. 
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Figure  51a.  Streamwise  turbulence  intensity  u2/l/|  profiles  at  X  =  4.27 
m,  k  =  0.61  flow  at  phases  120  0  ^  cot  ^  195°;  Legend  same  as  in  figure 
36b. 
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Figure  52A.  Streamwise  turbulence  intensity  profiles  at  X  =  4  27 

m,  k=1.33  flow  at  phases  195  0  £  cot  ^  255°;  Legend  same  as  in  figure 
39b. 
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Figure  53b.  Streamwise  turbulence  intensity  u2lU 2  profiles  at  X  =  3.59 
m,  k  =  1.03  flow  with  'roof  damper'  at  phases  105  °  <,  cot  <  180°;  Legend 
same  as  in  figure  40a. 
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Figure  53c  Streamwise  turbulence  intensity  u2/Ul  profiles  at  X  =  3.59 
m,  k  =  1.03  flow  with  'roof  damper7  at  phases  195  0  <,  a>t  <,  270°;  Legend 
same  as  in  figure  40b. 
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Figure  53d.  Streamwise  turbulence  intensity  u2lU2  profiles  at  X  =  3.59 
m,  k  =  1 .03  flow  with  'roof  damper'  at  phases  285  °  <,  cot  <  360°;  Legend 
same  as  in  figure  40b. 
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Figure  54a.  Normal-to-wali  turbulence  intensity  v2/Ul  premies  at  X 
4.27  m,  k  =  0.61  flow  at  phases  120  °  <,  wt  <,  195°;  Legend  same  as 
Figure  36b. 


Figure  55b.  Normal-to-wall  turbulence  intensity  profiles  at  X  = 

4.27  m,  k=1.33  flow  at  phases  270  °<  cof<  330°;  Legend  same  as  in 
figure  39b. 
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Figure  56a.  Normal-to-wall  turbulence  intensity  v2/U ]  profiles  at  X  : 
3.59  m,  k=1.03  flow  with  'roof  damper'  at  phases  15  °<  cot  <,  90 
Legend  same  as  in  figure  40a. 
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Figure  56b.  Normal-to-wall  turbulence  intensity  profiles  at  X  : 

3.59  m,  k  =  1.03  flow  with  'roof  damper'  at  phases  105  °  ^  (at  <,  180 
Legend  same  as  in  figure  40a. 
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Figure  56c.  Normal-to-wall  turbulence  intensity  v2/U|  profiles  at  X  =  3.59 
m,  k=  1.03  flow  with  'roof  damper'  at  phases  195  0  <,  cot  <,  270°;  Legend 
same  as  in  figure  40b. 
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Figure  56d.  Normal-to-wall  turbulence  intensity  v2/U2  profiles  at  X  : 
3.59  m,  k=1.03  flow  with  'roof  damper'  at  phases  285  0  ^  cot  <  360 
Legend  same  as  in  figure  40b. 
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e  57.  Normalized  Backflow  velocity  at  X  =  3.59  m,  k  =  1.03  flow  with 
damper';  Legend  same  as  in  figure  40a. 


y,  cm. 

Figure  60.  Variation  of  phase  angle  of  first  harmonics  at  X  =  2.95  m, 
k=  1.33  flow:  (  O  =  0:  A  =  -uv)\  open  symbols,  laser  anemometer; 
solid  symbols,  cross-wire  measurements. 


Figure  64a.-  uvmJU*  vsA/H  at  X  =  3.01  m„  k=1.03  flow  with  'roof 
damper',  at  X  =  3.59  m,  k  =  1 .03  flow  with  'roof  damper'.  Legend  same 
as  in  Figure  63a. 


Figure  64b.-  uvmjUl  vs.l/H  at  X  =  3.28  m.,  k=1.03  flow  with  'roof 
damper'.  Legend  same  as  in  figure  63a. 


Figure  67a.  \  UN\ /Ue  vs.  ]jH  for  k  =  0.61  flow,  solid  line  given  by  equation 
(11);  Dashed  Line  given  by  Equation  (12);  Leader  on  symbol  points  in 
the  local  direction  of  hysteresis  Loop  for  a  given  streamwise  location;  ( 
□  =  2.95 m;  +  =  3.51m;  A  =  4.27m  ). 
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Figure  67b.  \UN\/Uevs.  J\/H  for  k  =  1 .33  flow.  solid  line  given  by 
equation  (11);  Dashed  line  given  by  equation  (12);  Leader  on  symbol 
points  in  the  local  direction  of  hysteresis  Loop  for  a  given  streamwise 
location.  Legend  same  as  in  figure  67a. 


APPENDIX  A 


(la)  X 


cot 

(deg) 

15 

30 

45 

60 

75 

90 

105 

120 

135 

150 

165 

180 

195 

210 

225 

240 

255 

270 

285 

300 

315 

330 

345 

360 


Table  1  (k  =  0.61  flow) 

=  0.53  m 


A 

ue 

(m/s) 

A 

<5* 

(cm) 

A 

9 

(cm) 

4.84 

0.361 

0.236 

5.06 

0.312 

0.209 

6.62 

0.262 

0.182 

9.65 

0.204 

0.144 

13.24 

0.198 

0.153 

16.83 

0.195 

0.147 

20.50 

0.199 

0.155 

23.70 

0.221 

0.174 

25.60 

0.237 

0.190 

26.33 

0.224 

0.176 

26.75 

0.238 

0.186 

26.70 

0.220 

0.172 

n/-  nr\ 

0.228 

0477 

26.50 

0.232 

0.181 

26.24 

0.225 

0.175 

25.52 

0.200 

0.151 

23.60 

0.205 

0.158 

20.75 

0.219 

0.166 

17.50 

0.237 

0.178 

14.03 

0.260 

0.194 

9.74 

0.325 

0.229 

6.24 

0.454 

0.300 

4.62 

0.539 

0.346 

4.50 

0.465 

0.305 

Re a  H  Cfx  1 03 

722.7  1.530  5.356 

668.9  1 .494  5.203 

764.5  1.438  5.083 

876.9  1.417  5.085 

1280.3  1.297  4.656 

1567.9  1.325  4.424 

2009.0  1 .286  4.247 

2617.0  1-269  3.863 

3082.8  1.246  3.711 

2926.3  1.274  3.619 

3146.4  1.281  3.609 

2913.2  1.279  3.624 

2978.1  1.290  3.657 

3034.7  1.280  3.675 

2906.0  1 .284  3.746 

2443.3  1.322  3.798 

2358.6  1.301  3.828 

2185.7  1.315  3.918 

1976.8  1.331  3.995 

1726.9  1.336  4.121 

1409.1  1.421  4.160 

1184.2  1.513  4.225 

1010.7  1.55S  4. 

868.7  1.525  4. 


(deg) 

(m/s) 

15 

6.10 

30 

6.38 

45 

8.18 

60 

12.10 

75 

16.90 

90 

21.50 

105 

26.15 

120 

30.17 

135 

32.53 

150 

33.53 

165 

33.96 

180 

34.00 

195 

33.85 

210 

33.70 

225 

33.40 

240 

32.60 

255 

30.15 

270 

26.60 

285 

22.60 

300 

18.25 

315 

13.00 

330 

8.65 

345 

6.51 

360 

6.00 

1 1 9 

0 

75 

0 

51 

0 

H 

Cf  x  103 

4 

1.329 

4.371 

4 

1.331 

4.326 

6 

1.300 

4.428 

3 

1.255 

4.271 

4 

1.250 

3.783 

4 

1.242 

3.544 

7 

1.238 

3.318 

6 

1.240 

3.243 

4 

1.244 

3.160 

9 

1.241 

3.167 

5 

1.242 

3.224 

a 

1.246 

3.174 

9 

1.247 

3.180 

5 

1.246 

3.152 

3 

1.248 

3.250 

2 

1.251 

3.252 

3 

1.253 

3.341 

1 

1.249 

3.417 

2 

1.257 

3.490 

a 

1.278 

3.478 

4 

1.316 

3.460 

5 

1.352 

3.776 

* 

1.364 

3.868 

3 

1.343 

4.261 

V 
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1e)  X  =  2.95  m 


cot 

ue 

S' 

e 

Re a 

n 

(deg) 

(m/s) 

(cm) 

(cm) 

V 

90 

14.99 

2.408 

1.346 

12769.5 

105 

19.46 

2.229 

1.300 

16032.2 

120 

23.33 

2.321 

1.318 

19472.2 

135 

25.74 

2.436 

1.364 

22231.0 

150 

26.76 

2.492 

1.387 

23489.2 

165 

27.34 

2.547 

1.400 

24207.6 

180 

27.34 

2.517 

1.384 

23953.0 

195 

27.28 

2.544 

1.400 

24176.5 

210 

27.19 

2.485 

1.377 

23692.8 

225 

26.95 

2.456 

1.372 

23380.7 

240 

26.26 

2.542 

1.372 

22822.1 

255 

24.24 

2.724 

1.443 

22137.5 

270 

21.40 

3.049 

1.443 

19542.8 

285 

18.32 

4.189 

1.245 

14429.1 

14.90 

5.957 

0.278 

2618.0 

315 

10.65 

9.687 

-2.957  - 

19938.5 

330 

8.20 

13.989 

-7.117  - 

36943.5 

345 

7.56 

17.553 

-7.536  - 

36067.1 

1.789 

1.713 

1.761 


Cf  x  iO3 


0.909 

0.945 

0.812 


149 

(If)  X  =  3.51  m 


(Of 

(deg) 

A 

Ue 

(m/s) 

A 

<5* 

(cm) 

A 

0 

(cm) 

Re a 

0 

A 

H 

105 

15.16 

5.990 

3.321 

31869.7 

1.804 

120 

18.60 

6.160 

2.825 

33262.1 

2.181 

135 

20.76 

6.878 

2.816 

37005.0 

2.443 

150 

22.00 

7.274 

2.758 

38410.4 

2.638 

165 

22.62 

7.619 

2.675 

38304.7 

2.848 

180 

23.10 

6.961 

2.592 

37904.3 

2.686 

195 

23.11 

7.147 

2.603 

38078.7 

2.746 

210 

23.00 

6.946 

2.565 

37344.2 

2.708 

225 

22.85 

7.246 

2.545 

36809.1 

2.847 

240 

22.30 

7.583 

2.533 

35765.9 

2.993 

255 

20.82 

8.192 

2.511 

33102.3 

3.262 

270 

18.95 

9.288 

2.174 

26079.3 

4.273 

285 

16.88 

11.552 

1.752 

18717.5 

6.595 

300 

14.49 

15.911 

0.992 

9095.7 

16.047 

315 

11.30 

23.011 

-3.366 

-24079.1 

-6.836 

V 


J* 


ig)  x 

=  4.27  m 

A 

A 

A 

oot 

Ue 

<5* 

0 

Re a 

0 

A 

H 

(deg) 

(m/s) 

(cm) 

(cm) 

120 

12.70 

2.104 

5.465 

43939.3 

0.385 

135 

16.20 

8.370 

3.264 

33471.6 

2.565 

150 

18.50 

12.857 

2.710 

31739.5 

4.744 

165 

20.22 

17.105 

2.498 

31980.1 

6.847 

180 

21.10 

19.526 

2.625 

35060.0 

7.439 

195 

21.60 

19.778 

2.651 

36248.6 

7.461 

210 

21.90 

19.484 

3.307 

45855.3 

5.891 

225 

21.80 

19.542 

3.707 

51162.7 

5.271 

240 

21.25 

19.877 

3.607 

48524.1 

5.51  1 

255 

20.00 

20.640 

3.185 

40327.3 

6.480 

270 

18.40 

22.788 

2.320 

27020.1 

9.824 

28971.9 


2a)  X  =  0.53  m 

cot  Ue 
(deg)  (m/s) 


Table  2  (k  —  1.33  flow) 


A 

5' 


(cm) 


A 

e 


(cm) 


Re  a 

o 


A 

H 


VVVVV7/.J <V  «.  SJ.S  V 


2.95  m 

A 

A 

A 

ue 

<5‘ 

e 

Re a 

e 

m/s) 

(cm) 

(cm) 

7.80 

2.101 

1.271 

14317.5 

8.60 

2.285 

1.346 

15848.2 

8.70 

2.333 

1.294 

15319.3 

8.90 

2.551 

1.382 

16539.2 

8.63 

2.748 

1.442 

17001.5 

8.40 

2.819 

1.447 

16855.0 

8.11 

2.977 

1.491 

17092.0 

7.64 

3.159 

1.499 

16739.4 

7.08 

3.262 

1.520 

16437.3 

6.30 

3.437 

1.500 

15476.9 

4.65 

4.576 

1.308 

12130.9 

2.03 

6.296 

0.464 

3534.6 

9.48 

10.173 

-3.316 

-19901.0 

7.56 

18.442 

-14.842 

-71036.5 

6.45 

23.510 

-26.097 

-106563.1 

653 

0.942 

698 

0.827 

803 

0.887 

4 
9 

2.10 
2 


3.51  m 

A 

ue 

A 

<5* 

A 

0 

Re a 

e 

A 

H 

A 

Cf  x  ' 

(m/s) 

(cm) 

(cm) 

13.76 

4.071 

1.881 

16385.3 

2.164 

0.6i 

15.01 

4.850 

2.390 

22713.6 

2.029 

0.4! 

15.79 

6.401 

2.788 

27869.5 

2.296 

0.2< 

16.29 

7.470 

2.657 

27396.6 

2.812 

0.1! 

16.50 

8.798 

2.375 

24808.8 

3.705 

0.0: 

16.60 

10.498 

2.380 

25016.5 

4.410 

0.5: 

16.80 

11.585 

2.118 

22529.3 

5.469 

0.0 

16.64 

12.196 

1.815 

19118.2 

6.720 

o.o; 

16.10 

12.685 

1.662 

16943.5 

7.631 

0.1; 

19.06 

15.560 

1.904 

22968.7 

8.174 

O.t 

16.97 

16.194 

0.333 

3572.4 

48.699 

14.50 

19.475 

-1.922 

-17642.7 

-10.133 

11.30 

25.094 

-7.552 

-54021.8 

-3.323 

(2g)  X  =  4.27  m 


cot 

(deg) 

ue 

(m/s) 

(5* 

(cm) 

e 

(cm) 

Re« 

A 

H 

165 

10.12 

5.721 

3.540 

22678.4 

1.616 

180 

10.95 

5.186 

2.875 

19932.1 

1.804 

195 

11.91 

7.487 

3.360 

25332.6 

2.228 

210 

12.92 

11.052 

3.674 

30052.3 

225 

13.88 

15.744 

3.273 

28756.3 

4.811 

240 

14.60 

20.079 

2.119 

19582.3 

9.477 

255 

15.00 

23.954 

0.661 

6277.7 

36.235 

270 

15.17 

26.957 

-0.529 

-5079.3 

-50.968 

285 

14.50 

29.730 

-2.201 

-20206.4 

-13.506 

300 

12.84 

34.817 

-5.775 

-46945.0 

-6.029 

©SW5 


»J1*. 


159 


(3b)  X  =  3.01  m 


cot 

Ue 

<5‘ 

e 

CD 

H 

(deg) 

(m/s) 

(cm) 

(cm) 

15 

12.05 

5.748 

1.547 

11804.5 

3.714 

30 

12.30 

5.404 

1.661 

12933.8 

3.254 

45 

13.20 

4.373 

1.678 

14022.5 

2.606 

60 

14.30 

3.249 

1.456 

13184.7 

2.231 

75 

15.10 

2.695 

1.279 

12224.7 

2.108 

90 

15.50 

2.672 

1.217 

11942.4 

2.195 

105 

15.70 

2.837 

1.242 

12340.0 

2.285 

120 

15.90 

3.045 

1.258 

12659.4 

2.421 

135 

16.00 

3.210 

1.313 

13300.8 

2.444 

150 

16.00 

3.326 

1.341 

13586.6 

2.480 

165 

16.00 

3.386 

1.349 

13662.2 

2.510 

180 

16.00 

3.466 

1.381 

13991.8 

2.509 

195 

16.00 

3.410 

1.398 

14159.2 

2.440 

210 

16.00 

3.303 

1.340 

13574.5 

2.464 

225 

16.00 

3.285 

1.375 

13925.0 

2.390 

240 

15.90 

3.282 

1.352 

13608.6 

2.428 

255 

15.80 

3.371 

1.366 

13668.3 

2.467 

270 

15.70 

3.491 

1.369 

13607.8 

2.550 

285 

15.50 

3.622 

1.407 

13803.4 

2.575 

300 

15.28 

3.813 

1.428 

13814.0 

2.670 

315 

14.90 

3.872 

1.362 

12847.5 

2.843 

330 

14.30 

4.125 

1.356 

12274.7 

3.043 

345 

13.33 

4.499 

1.305 

11016.8 

3.446 

360 

12.30 

5.106 

1.279 

9955.9 

3.993 

ry; 


i!! 


. ..  . 
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(3c)  X  =  3.28  m 


cot 

Ue 

5' 

(deg) 

(m/s) 

(cm) 

15 

12.29 

12.452 

30 

12.53 

12.179 

45 

13.30 

1 1 .245 

60 

14.01 

9.621 

75 

14.36 

8.159 

90 

14.77 

7.863 

105 

15.26 

8.472 

120 

15.70 

9.172 

135 

15.92 

9.617 

150 

15.90 

9.689 

165 

15.90 

9.769 

180 

15.90 

9.775 

195 

15.91 

9.763 

210 

15.88 

9.618 

225 

15.80 

9.547 

240 

15.78 

9.679 

255 

15.70 

9.780 

270 

15.59 

10.206 

285 

15.41 

10.608 

300 

15.24 

10.855 

315 

14.95 

11.120 

330 

14.43 

11.396 

345 

13.60 

11.697 

360 

12.70 

12.222 

0  Re*  H 

u 


(cm) 

1.467 

11415.4 

8.487 

1.967 

15603.7 

6.191 

2.654 

22347.1 

4.237 

2.935 

26030.0 

3.278 

2.537 

23061.0 

3.216 

2.207 

20633.4 

3.563 

2.191 

21170.1 

3.866 

2.216 

22027.6 

4.138 

2.214 

22317.9 

4.343 

2.173 

21870.1 

4.459 

2.152 

21658.1 

4.540 

2.186 

22001 .0 

4.472 

2.223 

22389.5 

4.392 

2.183 

21943.9 

4.406 

2.220 

22208.8 

4.300 

2.251 

22482.4 

4.301 

2.140 

21270.0 

4.570 

2.111 

20838.1 

4.834 

2.010 

1 9606.8 

5.278 

1.983 

19133.5 

5.473 

1.954 

18492.4 

5.691 

1.863 

17014.6 

6.118 

1.648 

14187.9 

7.098 

1.386 

11142.5 

8.819 

■f  x  103 


3d)  X  =  3.59  m 


cot 

(deg) 

A 

Ue 

(m/s) 

A 

J* 

(cm) 

A 

0 

(cm) 

Re a 

6 

A 

H 

15 

12.27 

22.531 

-0.279 

-2166.4 

-80.787 

30 

12.38 

21.842 

0.487 

3817.2 

44.846 

45 

12.90 

22.021 

1.281 

10463.4 

17.187 

60 

13.37 

21.673 

2.103 

17802.2 

10.305 

75 

13.11 

19.791 

2.462 

20434.7 

8.038 

90 

13.14 

18.082 

1.615 

13438.5 

11.193 

105 

13.88 

17.970 

0.719 

6320.4 

24.983 

120 

14.87 

19.105 

0.450 

4236.3 

42.455 

135 

15.41 

20.228 

0./82 

7626.9 

25.874 

150 

15.37 

20.646 

0.868 

8445.4 

23.787 

165 

15.34 

20.607 

0.741 

7198.1 

27.802 

180 

15.43 

20.360 

0.573 

5595.6 

35.544 

195 

15.50 

20.493 

0.855 

8392.2 

23.961 

210 

15.43 

20.231 

0.894 

8732.5 

22.631 

225 

15.30 

20.208 

0.778 

7540.2 

25.958 

240 

15.20 

20.382 

0.520 

39.221 

255 

15.05 

20.975 

0.296 

2816.1 

70.966 

270 

14.85 

21.902 

0.037 

351.1 

586.413 

285 

14.64 

22.742 

-0.207 

-1921.3 

-109.705 

14.58 

23.532 

-0.154 

-1419.1 

-153.057 

315 

14.36 

23.545 

-0.162 

-1472.5 

-145.366 

330 

14.03 

23.515 

-0.212 

-1880.6 

-111.056 

345 

13.55 

23.279 

-0.183 

-1572.4 

-126.999 

360 

12.77 

23.134 

-0.388 

-3137.8 

-59.604 

aasa&as 
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APPENDIX  B 


DATA  ORGANIZATION 

Sets  of  the  ensemble-averaged  velocity  and  turbulence  quantities  data 
have  been  stored  on  diskettes  and  are  available  to  interested  users.  These 
data  are  stored  on  5.25  inch  high  capacity  floppy  diskettes,  which  can  be 
used  with  IBM/AT  with  high  capacity  drive  or  compatable  machines.  High 
capacity  diskettes  were  chosen  to  minimize  the  number  of  diskettes.  These 
data  may  become  available  also  on  60  mbyte  streaming  tape. 

Diskette  numbers  1  to  7  contain  data  for  k  =  0.61  and  1.33  flows  and  8 
to  12  contain  data  for  k  =  1.03  flow  with  "roof  damper";  as  shown  in  the 
following  table. 


Diskette 

number 

Shearwise 

location 

Flow 

Data  type 

(m) 

Code 

1 

0.53 

1 

k  =  0.61  &  1 .33 

U,u* , v2 ,-uv 

2 

1  .*15 

2 

ft 

tt 

3 

1.83 

3 

r» 

ft 

4 

2.20 

4 

tt 

ft 

5 

2.95 

5 

ft 

ft 

6 

3.51 

6 

ft 

tt 

7 

4.27 

7 

tt 

tt 

8 

1  .63 

1 

k  =  1.03  with 

"roof  damper" 

tt 

9 

3.01 

2 

tt 

tt 

10 

3.28 

3 

ft 

tt 

1 1 

3.59 

4 

ft 

ft 

12 

4.12 

5 

tt 

tt 

Data 

file  designation 

(i) 

First  letter  represents 

the 

flow. 

L  =  Low  frequency 

flow 

(f  =  0.596  Hz,  k  = 

0.61  ) 

H  <*  High  frequency 

flow 

(f  =  0.954  Hz,  k 

=  1  .33 

Total 
number  of 
data  files 


192 


96 


tt 


R  =  flow  with  "roof  damper"  (f  =  0.95*1  Hz,  k  =  1.03) 


163 


(ii)  Second  letter  represents  the  streamwise  distance  (see  the  table). 

(iii) Third  letter  (or  third  and  fourth)  letter  identify  the  type  of  data. 

U  =  ensembled-averaged  velocity 
U2  =  ensembled-averaged  u2 
V2  =  ensembled-averaged  v2 
UV  =  ensembled-averaged  -uv 

(iv)  Phase  angle  is  represented  by  last  two  (or  three,  as  the  case  may  be) 
digits. 

Example:  Data  file  H1U-180  represents  ensemble-averaged  velocity  data  for  k 
=  1.33  flow  at  X  =  0.53  m,  and  at  a  phase  angle  of  180  degrees. 

Data  file  organization 

All  data  files  contain  laser  and  valid  X-wire  data  except  for  R1(k  = 
1.03  flow  with  "roof  damper"  at  X  =  1.63  m),  where  only  X-wire  data  were 
taken.  First  line  of  the  data  files  contain  streamwise  location  (X  in 
inches),  phase  angle  (ut),  number  of  bins  (2M),  number  of  laser  points  (NL), 
number  of  X-wire  points  (NL)  and  a  dummy  number.  Next,  NL  line  contains 
laser  data  with  each  line  containing  four  entries,  i.e.,  distance  from  wall 

A  A  ^  ^  A  ^ 

in  Inches  (y),  U  (or  u2 ,  v2,  -uv),  uncertainty  ( AU  or  Au2  ,  av2,a-uv)  and  a 

~  A 

flag  (0,  for  laser  data).  X-wire  data  entries  contain  y,  U(or  u2,v2,-uv) 
and  a  flag  (1,  for  X-wire  data). 

A  program  to  read  the  data  files,  READ  is  available  on  each  diskette, 
READ  plots  the  laser  (small  circle)  and  X-wire  (larger  circle)  on  a  semi-log 
plot  on  the  screen.  X-axis  of  the  plot  is  the  distance  from  the  wall  on  a 
semi-log  scale  (3.6  cycle,  0.01  cm  to  60  cm)  and  Y-axis  is  the  quantity, 
i.e.,  U,u2,v2,  or  -uv.  Vertical  scale  is  scaled  between  minimum  and  maximum 
value  of  the  data  but  user  can  alter  the  scale  by  selecting  option  2.  User 


can  also  print  the  numerical  values  of  laser  and  X-wire  data  on  the  screen 
by  selecting  option  3  or  fl. 


